First-principles and semi-empirical calculations for bound hole polarons in KNDO3 
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The ab initio linear muffin-tin-orbital (LMTO) formalism and the semi-empirical method of the 
Intermediate Neglect of the Differential Overlap (INDO) based on the Hartree-Fock formalism are 
combined for the study of the hole polarons (a hole trapped nearby the cation vacancy) in a cubic 
phase of KNbC>3 perovskite crystals. The 40-atom and 320-atom supercells were used, respectively. 
We predict existence of both, one-site and two-site (molecular) polarons with close optical absorption 
energies (0.9 eV and 0.95 eV). The relevant experimental data are discussed. 
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I. INTRODUCTION 

It is well understood now that point defects play an 
important role in the electro-optical and non-linear op- 
tical applications of KNLO3 and related ferroelectric 
materials.El In particular, reduced KNb03 crystals con- 
taining oxygen vacancies reveal fast photorefractive re- 
sponse to short-pulse excitations which could be used for 
developing fast optical correlators.0 The prospects of the 
use of KNb03 for the light frequency doubling are seri- 
ously affected by the presence of unidentified defects re- 
sponsible for induced infrared absorption.El The photore- 
fractive effect, important in particular for holographic 
storage, is also well known to depend on the presence 
of impurities and defects. Most of as-grown ABO3 per- 
ovskite crystals are non-stoichiometric and thus contain 
considerable amount of vacancies. 

The so-called F and F + centers (an O vacancy-, V , 
which traps two or one electron, respectively) pu be- 
long to the most common defects in oxide crystals. In 
electron-irradiated KNLO3, a broad absorption band ob- 
served around 2.7 eV at room temperature has been ten- 
tatively ascribed to the F-type centersE These two de- 
fects were the subjectjaf our recent ab initio and semi- 
empirical calculations .Em A transient optical, absorption 
band at 1.2 eV has been associated recently)I3 in analogy 
with other perovskites, with a hole polaron (a hole bound, 
probably, to a K vacancy). The ESR study of KNb0 3 
doped with Ti 4+ gives a proof that holes could be trapped 
by such negatively charged defects.til For example, in 
BaTiC>3, the hole polarons bound to Na and K alkali ions 
replacing Ba andr-fchus forming a negatively charged site 
attracting a holeJl3 have also been found. Cation vacan- 
cies are the most likely candidates for pinning polarons. 
In irradiated MgO, they are known to trac-one or two 
holes giving rise to the V~ and V° centerafrQ which are 
nothing but bound hole polaron and bipolaron, respec- 
tively. The results of the experimental studies of hole 
polarons in alkali halides and ferroelectric perovskites re- 
veal two different forms of atomic structure of polarons: 



atomic one (one-site), when a hole is localized on a sin- 
gle atom, and molecular-type (two-site), when a-hple is 
shared by two atoms forming a quasi-molecule.EnI3 In 
the present study, we simulate hole polarons associated 
with a K vacancy in KNb03, using an ab initio density 
functional theory (DFT)-based method in combination 
with a semi-empirical treatment based on the Hartree- 
Fock (HF) formalism, employing the periodic boundary 
conditions and the supercell geometry in both cases. 



II. METHODS 

The motivation for using the DFT-based and HF-based 
calculation methods in parallel is to combine strong sides 
of both in a single study. The DFT is expected to be able 
to provide good description of the ground state, i.e. to 
deliver reasonable relaxation energies and ground-state 
geometry. In the HF approach, the relaxation energies 
are generally less accurate because of the neglection of 
correlation effects. On the other hand, the HF formalism 
is straightforwardly suited for the evaluation of excitation 
energies, because the total energies can be calculated for 
any (ground-state or excited) electronic configuration on 
equal footing, that is generally not the case in the DFT. 
Practical experience shows that HF and DFT results of- 
ten exhibit similar qualitative trends in the description 
of dielectric properties but quantitatively lie on opposite 
sides of experimental data, thus.-effectively setting error 
bars for a theoretical prediction.^ 

Our ab initio DFT treatment is based on the full- 
potential linearized muffin-tin orbital (LMTO) formal- 
ism, previously applied with success to the study 
of structural instability and lattice dynamics in pure 
KNb03.E3il3 For the study of defects, we used the ver- 
sion of LMTO as implemented by van Schilfgaarde and 
McthfcssclEJ, that was earlier used in our simulations 
for the F-center in KNb03 (Ref. ^ see also for more 
details of the calculation setup there). The exchange- 
correlation has been treated in the local density approxi- 
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mation (LDA), as parametrized by Perdew and Zunger.La 
The supercell LMTO approach has been earlier success- 
fully j^ed for the simulation of defects in KCllla and 
MgO.tHl In the present case we used the 2x2x2 super- 
cells, i.e. the distance between repeated point defects was 
pb8 A. As a consequence of the large number of eigen- 
states per k-point in a reduced Brillouin zone (BZ) of 
the supercell and of the metallicity of the doped system, 
it was essential to maintain a quite dense mesh for the 
k-integration by the tetrahedron method over the BZ. 
Specifically, clear trends in the total energy as function of 
atomic displacements were only established at 10x10x10 
divisions of the BZ (i.e., 186 irreducible k-points for a 
one-site polaron). 

In the HF formalism, a semi-empirical I 
Neglect of the Differential Overlap (IND. 
modified for ionic and partly ionic solids ,c3 
used. The INDO method has been success 
for the simulation of defects in many oxides.l 
calculations have been performed with periodic boundagz. 
conditions in the so-called large unit cell (LUC) model,E3 
i.e., for k=0 in the appropriately reduced BZ. When 
mapped on the conventional BZ of the compound in ques- 
tion, this accounts for the band dispersion effects and al- 
lows to incorporate effectively the k-summation on a rela- 
tively fine mesh over the BZ. Due to the robustness of the 
summation procedure, the total energy dependence on 
the atomic displacements was found satisfactorily smooth 
for any supercell size (one should keep in mind, however, 
that the numerical results may remain somehow depen- 
dent on the LUC in question). In the present work, 4x4x4 
supercells (320 atoms) were used, that is the extension 
of our preliminary INDO studyO with the 2x2x2 super- 
cells. The detailed analysis of the INDO parametrisation 
for KNb03 is presented in Ref. |3(], and the application 
of the method to the F-center calculations - in Ref. |]|]. 

We restricted ourselves to a cubic phase of KNDO3, 
with the lattice constant ao=4.016 A. In the vacancy- 
containing supercell, the relaxation of either one (for the 
one-site polaron) or two neighboring (for the two-site po- 
laron) O atoms, amongst twelve closest to the K vacancy, 
has been allowed for, and the changes in the total energy 
(as compared to the unrelaxed perovskite structure with 
a K atom removed) have been analyzed. Also, we stud- 
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FIG. 1. Uniform relaxation of O atoms around the K va- 
cancy in KNDO3 (a) and the relaxation corresponding to the 
formation of one-site (b) and two-site (c) polarons. Arrows in 
(b) and (c) indicate the displaced O atoms. 
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FIG. 2. Local DOS at the K vacancy site and at the adja- 
cent oxygen atom (top panel) and at Nb and O sites in perfect 
KNb03 (bottom panel), as calculated by LMTO. 



ied the fully symmetric relaxation pattern (breathing of 
twelve O atoms) around the vacancy. Different relax- 
ation patterns considered in the present study are shown 
in Fig. [l]. The positions of more distant atoms in the 
supercell were kept fixed. 



III. RESULTS 

The removal from the supercell of a K atom with its 7 
electrons contributing to the valence band (VB) produces 
slightly different effects on the electronic structure, as de- 
scribed within the DFT and in the HF formalism. Acord- 
ing to the LMTO result, the Fermi energy lowers, and 
the system becomes metallic (remaining non- magnetic). 
Therefore, no specific occupied localized state is associ- 
ated with the vacancy. The local density of states (DOS) 
at the sites of interest is shown in Fig. |[ As is typi- 
cal for LDA calculations, the one-electron band gap in 
KNb03 comes out underestimated («2 eV) as compared 
to the experimental optical gap («3.3 eV). The removal 
of a K4s electron amounts to adding a hole which forms 
a localized state at «10 eV above the Fermi level, i.e. 
above the unoccupied Nb4d band. In the 2p-DOS of O 
atoms neighboring the vacancy a quasi-local state (that 
effectively screens the hole) is visible just below the Fermi 
level. Apart from that, the 02p-DOS is largely unaffected 
by the presence of vacancy and the changes in the DOS 
of more distant sites (K, Nb) are negligible as compared 
with those in the perfect crystal. As the cubic symmetry 
is lifted by allowing a non-uniform relaxation of O atoms, 
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the "screening" quasi-local state is clearly localized at the 
atom closest to the vacancy. At the same time, the hole 
state becomes smeared out in energy. This amounts to 
the bonding being established between the hole and the 
screening charge on one of its neighbors. 

In the INDO treatment, the one-electron optical gap 
is overestimated, as is typical for the HF calculations 
(«6 eV, see Ref. but the ASCF gap for the triplet 
state is 2.9 eV, close to the experiment. The quasi-local 
"screening" state is described by a wide band close to 
the VB top. This is consistent with the LDA description. 
But the INDO calculation also suggests, and this differs 
from the LDA, that the removal of an electron leaves 
an unpaired electron state split-off at ssleV above the 
VB band top. In case of asymmetrical O relaxation, the 
molecular orbital associated with this state is centered at 
the displaced O atom, being a combination of the 2p x , 
2p y functions of the latter in the setting when the plane 
spanned by their lobes passes through the vacancy site. 
The same applies qualitatively to the two-site polaron, 
with the only difference that the localized state is formed 
from the 2p orbitals of both O atoms approaching the 
vacancy, with a corresponding symmetry lowering. The 
localized hole state is also present in the HF description 
but lies much lower than the corresponding state in the 
LDA, forming a 0.9 eV -wide band located w 0.2 eV below 
the conduction band bottom (see Fig. |3j) . 

Differently from the DFT-based approaches which ad- 
dress in principle only the ground-state electron density, 
the HF method provides a possibility to evaluate exci- 
tation energies by means of the so-called ASCF formal- 
ism, i.e. as the difference of total energies in relaxed 
ground-state and excited states. In agreement with the 
Schirme|r's theory for the small-radius polarons in ionic 
solids,Eil the optical absorption corresponds to a hole 
transfer to the state dclocalizcd over nearest oxygens. 
The absorption energies due to the electron transition 
from the quasi-local states near the VB top (1, Fig. |^) 
into the vacant polaron band (2, Fig. ||) for one-site and 
two-site polarons are close (Table I), and both are twice 
smaller than the experimental value for a hole polaron 
trapped by the Ti impurity.liil This shows that the op- 
tical absorption energy of small bound polarons can be 
strongly dependent on the defect involved. Another im- 
portant observation is that the ASCF energy for absorp- 
tion turns out to be considerably smaller than the esti- 
mate based on the difference of one-electron energies. 

In spite of generally observed considerable degree of 




FIG. 3. Sketch of the polaron optical transition from the 
quasi-local state 1 near the top of the valence band to the 
hole state 2 below the conduction band bottom. 3 indicates 
the level of an unpaired electron. 
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FIG. 4. Relaxation energy for one-site polaron according 
to LMTO and INDO calculations vs the displacement of one 
oxygen atom towards the vacancy (a); relaxation energy for 
two-site polaron as function of the center-of-mass displace- 
ment of the O-O pair and of the O-O distance according to 
the INDO calculation (b). 

covalency in KNb03 and contrary to a delocalized char- 
acter of the F center statejm the one-site polaron state 
remains well localized at the displaced O atom, with only 
a small contribution from atomic orbitals of other O ions 
but none from K or Nb ions. Although there are some 
differences in the description of the (one-particle) elec- 
tronic structure within the DFT- and HF-based meth- 
ods, the trends in the total energy driving the struc- 
ture optimization remain essentially the same. In both 
approaches, both one-site and two-site configurations of 
the hole polaron are much more energetically favorable 
than the fully symmetric (breathing mode) relaxation of 
twelve O atoms around the K vacancy. This is in line 
with what is,ktiown about small-radius polarons in other 
ionic solidsEjtiJ and is caused by the fact that the lattice 
polarization induced by a point charge is much larger 
than that due to a delocalized charge. 

In the case of one-site polaron, a single O - ion is dis- 
placed towards the K vacancy by 1.5 % of the lattice con- 
stant (LMTO) or by 3% (INDO) - see Fig. |. The INDO 
calculations show that simultaneously, 11 other nearest 
oxygens surrounding the vacancy tend to be slightly dis- 
placed outwards the vacancy. In the two-site (molecular) 
configuration, a hole is shared by the two O atoms which 
approach each other - by 0.5% (LMTO) or 3.5% (INDO) 
- and both shift towards a vacancy - by 1.1% (LMTO) 
or 2.5% (INDO). The lattice relaxation energies (which 
could be associated with the experimentally measurable 
hole thermal ionization energies) are presented in Table 



TABLE I. Absorption (E^^) and lattice relaxation (£" re p 
relatively to the perfect crystal with a K vacancy) energies 
(eV), calculated by LMTO and INDO methods. 



Energy: 


E a i, B 


E re i 


E re ; 


Method: 


INDO 


LMTO 


INDO 


uniform breathing 




0.01 


0.08 


one-site polaron 


0.90 


0.12 


0.40 


two-site polaron 


0.95 


0.18 


0.53 
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I. In both methods the two-site configuration of a polaron 
is lower in energy. 

A comparison of the present, 320-atomic INDO calcu- 
lation with a preliminary calculationE3 using 40— atomic 
LUC and self-consistency in the T point of the BZ only 
shows that the optical absorption energies are changed in- 
considerably, unlike the lattice relaxation energies. The 
latter now are much smaller and thus in better agreement 
with the LMTO calculation. 



IV. CONCLUSIONS 

In this pilot study we focused on the quantitative mod- 
els of hole polarons in KNbC>3. The main conclusion is 
that both one-center and two-center configurations are 
energetically favorable and close in energy (with a slight 
prefernce for the two-center configuration), as follows 
from the numerical simulation results by two different 
theoretical methods. The calculated optical absorption 
energies and the spatial distribution of relevant electronic 
states could provide guidelines for more direct exper- 
imental identification of the defects in question. The 
calculated hole polaron absorption («1 eV) iSj-close to 
the observed short-lived absorption band energyE3; hence 
this band could indeed arise due to a hole polaron bound 
to a cation vacancy. Further detailed study is needed 
to clarify whether such hole polarons are responsible 
for the effect of the blue light-induced infrared absorp- 
tion reducing the second-harmonic generation efficiency 
in KNb0 3 Ja 

As compared with the DFT results, the INDO (as is 
generally typical for the HF-based methods) systemat- 
ically gives larger atomic displacements and relaxation 
energies. The DFT results are more reliable in what re- 
gards the ground state of polarons, whereas the use of 
the HF formalism was crucial for the calculation of their 
optical absorption and hence possible experimental iden- 
tification. 
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